How periodic orbit bifurcations drive multiphoton ionization 
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The multiphoton ionization of hydrogen by a strong Dichromatic microwave field is a complex 
process prototypical for atomic control research. Periodic orbit analysis captures this complexity: 
Through the stability of periodic orbits we can match qualitatively the variation of experimental 
ionization rates with a control parameter, the relative phase between the two modes of the field. 
Moreover, an empirical formula reproduces quantum simulations to a high degree of accuracy. This 
quantitative agreement shows how short periodic orbits organize the dynamics in multiphoton ion- 
ization. 
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Simple systems can display extraordinarily complex 
dynamics: This lesson from three decades of chaos the- 
ory has altered the direction of many areas of physics 
One-electron systems, which are the most funda- 
mental at the atomic level, have had their fair share of 
such striking new discoveries Q, a prominent example 
being the multiphoton ionization of hydrogen in a strong 
microwave field [|[. Its interpretation remained a puz- 
zle until its stochastic, diffusional nature was uncovered 
through the then- new theory of chaos Feverish re- 
search activity in the last thirty or so years has resulted in 
a fairly complete understanding of this problem [5J, [y, [7| • 
As in many other areas of physics, here too the empha- 
sis has shifted lately from understanding the process to 
using those insights to manipulating it [1, Q. Clearly, 
the hope is that experience gained from this prototyp- 
ical system will help to control more complex systems 
ranging from atoms to plasmas. Here control denotes 
tailoring the physical behavior of nonlinear dynamical 
systems (which generically exhibit chaotic dynamics) us- 
ing "knobs" (i.e., suitable external parameters). Success 
depends on identifying simple knobs and understanding 
why and how they alter the system. In this Letter we 
report on precisely such a knob for a complex quantum 
system, and show how ionization behavior can be pre- 
dicted with quantitative accuracy using periodic orbits. 

The ionization of a one-dimensional hydrogen atom 
driven by a bichromatic linearly polarized electric field 
is a seemingly simple set-up with complex dynamics. 
Bichromatic pulses [1(| HH have emerged as natural tools 
in atomic control research because the relative phase 
of the two phase-locked pulses offers a practical control 

knob m m m ee us m m m m m m. The 

Hamiltonian of the system is, in atomic units, 
P 2 1 

H = h F h x sinChut) + Fixsmtltot + 0), (1) 
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where the indices I , h refer to the low and high frequency 
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modes with frequencies Ilj and hu;, respectively. These 
two modes are frequency locked, i.e. I and h are integers. 
They are out-of-phase by <f), the control parameter. 

In what follows, we consider Hamiltonian (Q} for the 
two sets of experimental parameters of Ref . [20| : (/) the 
h:l=3:l, F h = 24 VcuT 1 and Fi = 53.4 Vcm" 1 , and (77) 
the h:l=3:2, F h = 25 Vcm" 1 and F t = 33.5 Vcm" 1 . The 
highest frequency is 18 GHz in both cases. These two 
sets show experimentally and numerically drastically dif- 
ferent ionization behavior [2(| and in the frequency range 
considered, these results can be reproduced by quantum 
or classical simulations. Our purpose here is to show 
that these findings can be qualitatively and quantita- 
tively captured using a periodic orbit analysis, which re- 
veals the classical bifurcations responsible for ionization. 
It also allows prediction of ionization at other values of 
parameters without resorting to large numerical simula- 
tions. 

We begin by mapping Hamiltonian ([1]) into action- 
angle variables such that the principal quantum number 
n is associated with action J. We assume uj = 1 without 
loss of generality. Hamiltonian ([I]) becomes Q 

H = -^ + 2J 2 [F h sm(ht) 

+F t ahx(lt + <j))] + Y2 ak cos k9 ^j > ( 2 ) 

where a n = [J n (n) — J n _i(n)]/n and J„'s are Bessel func- 
tions of the first kind. Note that there are three variables 
(J, 0, t) and three parameters (Fh,Fi,<j>). We denote this 
Hamiltonian H(J,0,t;4>) where (Fh,Fi) are chosen ac- 
cording to Cases (I) or (//). 

For this Hamiltonian, we consider a specific periodic 
orbit, denoted O(0) for <fi — 0. Numerically it is deter- 
mined using a modified Newton-Raphson multi-shooting 
algorithm as described in Ref. [l| . We follow this periodic 
orbit as cj> is varied. The orbit O(0) deforms continuously 
into 0(4>), the period of which is denoted T(<f>). In ad- 
dition to its location, we also monitor its linear stability 
properties obtained from the reduced tangent flow ex- 
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pressed as 

where J is the two-dimensional skew-symmetric matrix 
and V 2 £f is the two-dimensional Hessian matrix (com- 
posed by second derivatives of H with respect to its 
canonical variables J and 6). The initial condition is 
J® = I2 (the two-dimensional identity matrix). The 
stability properties are given by the two eigenvalues of 
the monodromy matrix jl which form a pair (A, 1/A) 
(since the flow is volume preserving, the determinant of 
is equal to 1). The periodic orbit is elliptic if 

the spectrum is (e 1 "^, e~ w ^) (and stable, except at 
some particular values), or hyperbolic if the spectrum is 
(A((j6), 1/X((f>)) with A(0) € R* (unstable). In a more con- 
cise form, it can be summarized using Greene's residue, 

r 0, m 

^) = (2-trjf«)/4. 

If R{4>) G]0, 1[, the periodic orbit is elliptic; if R(<fi) < or 
R(4>) > 1 it is hyperbolic; and if R(4>) = and R{4>) = 1, 
it is parabolic. Generically, periodic orbits and their lin- 
ear stabilities are robust against small changes of param- 
eters, except at specific values where bifurcations occur 
[25l ]. These rare events affect the dynamical behavior 
drastically. 

In what follows, we identify the bifurcations (if any) of 
a set of short periodic orbits, i.e., the type and the value 
of the parameter cj) c where they bifurcate. This provides 
a way to foretell if a relatively high ionization rate should 
be expected or not. We use the residue curves 4> l— ► -R(</>) 
for each periodic orbit to analyze the dependence of ion- 
ization rates on <p . The importance of considering two 
associated Birkhoff periodic orbits (i.e. periodic orbits 
with the same action but different angles in the integrable 
case, one elliptic and one hyperbolic), was emphasized in 
Ref. [H 

Figure Q] shows a Poincare section of Hamiltonian @ 
for Case (I) at <fi — 0. We notice two main islands in 
the chaotic sea. At the centers of these islands sit elliptic 
periodic orbits with period 27r (indicated by full circles). 
Note that the (rescaled) principal quantum number con- 
sidered in Ref. [2(| lies in between these two islands. The 
residue method will monitor these two elliptic periodic 
orbits as well as their associated hyperbolic periodic or- 
bits (they all have the same period T(<fi) = 2ir for all 
values of (f>) [26| . We follow the location and residue for 
each of these periodic orbits as the parameter <j) varies. 
Figure [5] shows the residue curves of these orbits. 

In Fig. [5] (upper panel) , if we start from </> = and fol- 
low the upper elliptic periodic orbit of Fig. [TJ it remains 
elliptic (R(4>) G]0, 1[) until </> c ss 0.49 where a bifurcation 
appears. At this critical value, the orbit is parabolic. 
Then it turns hyperbolic {R{4>) > 1) until 27r/3— 4> c where 
another bifurcation occurs and the orbit returns to being 




FIG. 1: Poincare section of Hamiltonian ([2]) for Case (J) at 
4> — 0. Full circles and crosses indicate the two elliptic and 
two hyperbolic periodic orbits we consider respectively. The 
horizontal line corresponds to the principal quantum number 
n = 51 considered in Ref. [2pj | . 
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FIG. 2: Residue curves for the four periodic orbits with period 
2tt indicated by crosses and circles on Fig. [l] The bold curves 
are for the upper set of elliptic/hyperbolic orbits. Small ar- 
rows indicate where the bifurcations happen. Upper panel is 
for Case (7), i.e. h:l—3:l and lower panel for Case (//), i.e. 
h:l=3:2. 

elliptic. This bifurcation is of the period doubling kind 
at (j) c and a period halving one at 27r/3 — cf> c . The compu- 
tation of log|A±(0)| (where \±((f>) are the two eigenval- 
ues of the monodromy matrix J'T^ associated with the 
upper elliptic periodic orbit) show that, before the bifur- 
cation, log|A(0)| = since the orbit is elliptic, and that 
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just after the bifurcation, log |A±(0)| oc — (/> c - Note 
that while the upper elliptic periodic orbit undergoes a 
bifurcation as <fi is varied, the other three periodic orbits 
retain the stability properties they had at (j> = 0. 

In the parameter range <f> £]</> c ,27r/3 — (f> c [, the up- 
per part of phase space does exhibit more chaos. Since 
the initial atomic beam is taken in this region (principal 
quantum number n — 51 which corresponds to action 
Ji s» 0.495), the ionization rate is expected to be higher 
in this regime. Apart from fine detail (like higher order 
islands), the upper part of phase space is roughly homo- 
geneous. Therefore we expect to observe a plateau in 
the ionization probability versus <p. The reason is that 
a strongly hyperbolic orbit only influences the ionization 
time and not the value of the ionization probability. Of 
course, this is true provided that the duration of the max- 
imum pulse envelope is large enough (in the experiment, 
this is approximately 15 times the period of the periodic 
orbits considered |20||). Roughly speaking, this means 
that in the chaotic region all the orbits ionize (i.e., es- 
cape to a value of the action J lon > 1.26) whatever the 
hyperbolicity degree is. In Ref. [20], experimental results 
as well as one-dimensional quantum calculations show 
this plateau. Fro m q uantum calculations, </> c rs 0.5 was 
obtained in Ref. [20| which is in good agreement with 
the parameter value <pc ~ 0.49 at which the bifurcation 
of the upper elliptic periodic orbit occurs. 

We perform the same analysis for Case (77) where 
h:l=3:2. Again we consider two main islands in the 
chaotic sea where two different elliptic periodic orbits 
with period 2-7T sit at the centers. We monitor the stabil- 
ity of four periodic orbits (two elliptic orbits and the as- 
sociated hyperbolic ones) . Figure [2] (lower panel) shows 
the four residue curves. Evidently the elliptic periodic 
orbits remain elliptic and the hyperbolic ones remain hy- 
perbolic for all values of 0. No bifurcations occur; conse- 
quently, the ionization probability is expected to be ap- 
proximately independent of and to be lower than Case 
(7) since for these values of amplitudes, the chaotic re- 
gion is smaller. This is consistent with the experimental 
and quantum calculations of Ref. [201 ( see their Fig. 7). 
The experimental results show a nearly flat curve for the 
ionization probability versus <p, whereas the quantum cal- 
culations show significant variations for this probability 
but no sharp increase and decrease as in Case (I). 

The periodic orbit analysis above elucidates whether or 
not there is a significant ionization probability for specific 
parameter values, and also where plateaus are expected 
to occur. This qualitative agreement highlights the im- 
portant role played by these orbits. Furthermore, we 
can obtain quantitative agreement concerning the shape 
of the ionization curve versus the phase <fi by using the 
residue curves. To this end, we devise an empirical for- 
mula for relative ionization probability in the following 
way : First, the values of 4> giving' the highest ionization 
would be the ones associated with the highest variations 
of the residues (in absolute value) with respect to the 
minimum ionization. Second, if the periodic orbit is too 
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FIG. 3: Normalized ionization rate vs (f> based on Eq. ((3]) for 
Case (I) (with fc:Z=3:l) with A = -2.52 and B = 0.65 (upper 
panel), and for Case (II) (with h:l=3:2) with A = -0.485 
and B = 0.17. Circles represent the data obtained by one- 
dimensional quantum calculations, taken from Ref. 
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far (in action) from the considered action J; then it will 
not influence the dynamics so there should be a penal- 
izing term depending on its position with respect to the 
chosen rescaled action. The formula reads : 



-Pion(0) =A + B 



M 



exp|7? m (0) -R m (0)\ 
exp \J m (4>) - Ji\ 



(3) 



where the sum is taken over the M different periodic 
orbits considered and J m (</>) = J(6)d6/2 / n is the av- 
erage action of the periodic orbit m. The parameters A 
and B in Eq. ([3]) are merely a translation and a dilata- 
tion of the curve in order to match the mean value and 
the amplitude of variations of P mri obtained in Ref. [2(| • 
This formula takes into account the value of the residues 
at <j) = 0. The aim is to set up a baseline for each of the 
periodic orbits (which is taken here at the value of the 
parameter where the ionization is minimal [20l]). In gen- 
eral, Eq. ([3]) can exhibit values which are greater than 
1, which are not relevant. In order to remedy to this 
problem, we truncate P; on at the value where a bifur- 
cation occurs. Therefore in the range where R n (4>) is 
larger than one, Pi on is constant (taken as the value of 
the residue at (j) c where the bifurcation occurs). 

Figure[3]depicts P; on given by Eq. ([3]) versus parameter 
4> as well as the data taken from Ref. 20] for both cases 
[Case (I) (upper panel) and Case (77) (lower panel)]. 
Since there are no bifurcations in Case (77), there is no 
plateau. We notice that the empirical formula reproduces 
accurately the results obtained from quantum calcula- 
tions. In particular, it captures some essential features 
of the ionization curve in Case (77) , like the two unequal- 
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FIG. 4: Bifurcation surface in parameter space (cj>,Fh,Fi) for 
h:l=3:l. 



TABLE I: Ionization thresholds obtained for Fh — 6 Vcm - , 
experimentally in Ref. pj| and by the residue method (see 
Fig- HI- The If case corresponds to Fh = 0. 



Fi{ Vcm- 1 ) 


cf>= 


4> = n/3 


If 


[18] 


107 


85 


96 


residue 


109.6 


81.4 


96.2 



sized peaks and the specific shape of both peaks. 

We can also predict the behavior of the system as all 
three parameters (the two amplitudes and the relative 
phase of the field) are varied. In Fig. [4] we represent the 
set of parameters where the upper elliptic periodic orbit 
of Fig. Q] (in the 3:1 case) is in fact parabolic (i.e., the set 
of parameters where the system undergoes a major bifur- 
cation). The equation of this surface in parameter space 
is R{4>, Fh,Fi) — 1. The boundaries of the plateaus in pa- 
rameter <p of Fig. [3] (upper panel) obtained by fixing the 
two values for Fh and Fi are on this surface. We notice 
that as expected, when Fh approaches zero, this surface 
is less dependent on parameter <f>. Table U com par es the 
experimental values for ionization thresholds [18l ] with 
corresponding values from our bifurcation analysis. 

To obtain this ionization surface in parameter space, a 
large number of classical trajectories for each value of the 
parameters (4>,Fh,Fi) needs to be computed for a suffi- 
ciently long time in order to decide whether or not a given 
trajectory has ionized. In contrast, only one orbit for a 
short time (typically the period of the field) is needed for 
the residue analysis. Furthermore, using residues, this 
surface can be constructed locally without any need to 
consider all possible values of the parameters. 

More broadly, our results constitute rules by which 
this quantum system can be controlled. Such systematic 
and practical coherent control rules remain very rare and 
sought after 
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